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Abstract 

We study the impacts of the satelhte galaxies in a measurement of the redshift- 
space distortions. In our multipole power spectrum analysis of the luminous red 
galaxies (LRGs) samples in the Sloan digital sky survey (SDSS), we have clearly de- 
tected the non-zero signature of the hexadecapole and tetrahexadecapole spectrum, 
which disappears in the power spectrum with the sample with the brightest LRGs 
only. We thus demonstrate that the satellite LRGs in multiple systems significantly 
contribute to the multipole power spectrum though its fraction is small. The behav- 
ior can be understood by a simple halo model, in which the one-halo term, describing 
the Finger of God effect from the satellite galaxies, is the main source of the higher 
multipole spectra. We also discuss about a prospect for the contribution of satellite 
galaxies in future large redshift surveys targeting the H-alpha emitters. 



1 Introduction 



The luminous red galaxies (LRGs) in the Sloan digital sky survey (SDSS) demonstrated 
the usefulness of a large redshift-survey of galaxies. Especially, it proved that a precise 
measurement of the statistical features in their spatial distribution provides us with the 
very useful methodology not only for the cosmology but also for the fundamental physics. 
For example, the baryon acoustic oscillation signature in the large scale structure is now 
recognized as a promising way for exploring the origin of the accelerated expansion of 
the universe [1, 2]. A stringent constraint on the neutrino mass is also obtained [3, 4]. 
Furthermore, the LRG sample showed that a measurement of the redshift-space distortions 
gives us a unique chance of testing the theory of gravity (e.g., [5, 6], cf. [7]). 

The LRGs in SDSS are massive early-type galaxies, and most part of them are consid- 
ered to be residing in the center of halos. However, it is clarified that the some fraction of 
the LRGs consists of multiple galaxies system. The halo occupation distribution (HOD) 
of the LRGs was clarified by Reid and Spergel [8]. The correspondence between the LRGs 
and halos has been investigated, which has illuminated the importance of the Finger of 
God (FoG) effect [3, 9, 10, 11]. Recent investigations with N-body simulations have dis- 
covered that halos at the redshift 2 could be the origin of the LRG host halos [12]. In the 
present paper, we investigate the contribution of the satellite LRGs in multiple systems 
to the redshift-space distortions. A related topic has been investigated in the literature 
[3, 9, 10, 11], but the previous works investigated the contribution to the monopole spec- 
trum. We here focus our investigation on the redshift-space distortions described by the 
higher multipole power spectrum. 

The redshift-space distortions are measured in terms of the anisotropic correlation 
function or the anisotropic power spectrum, e.g., [13, 14, 7]. The anisotropic correlation 
function of the SDSS LRG sample has been measured in the literature, e.g., [2, 15, 16]. 
The anisotropic power spectrum P{k,fi), where fj, denotes the directional cosine between 
the line of sight direction and the wave number vector, is the Fourier transform of the 
anisotropic correlation function. They are equivalent to each other. The multipole power 
spectrum P£{k) is defined as the coefficient of the multipole expansion,*^ *^ 



In refs. [5, 19], the multipole power spectrum of the SDSS LRG sample was measured. In 

"^^Note that the definition of Pe{k) is different from the conventional definition by tlie factor 2£ + 1. 
"^^There are also different works using phase-space distribution function approach to study the redshift 
distortion effect on the multipole power spectrum [17, 18]. 



P{k,fi) = Y,Pi{k)W{2i+l), 01 P,{k) = - P{k,fi)C,{^i)dfi, (1) 



where Ce{fi) is the Legendre polynomial, which is normalized as 




(2) 
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the present paper, we demonstrate that the satelhte galaxies significantly contribute to 
the higher multipole power spectrum, though its fraction is small. 

The primary purpose of the present paper is to understand the contribution of satellite 
galaxies to the multipole power spectrum. To this end, we measure the multipole power 
spectrum of the SDSS LRG samples, and compare the results with the predictions of a 
simple halo model with the HOD of the SDSS LRG catalog. Then, we show the importance 
of the one halo term in the higher multipole power spectrum. We also investigate the 
influence of satellite galaxies in a future redshift survey on the multipole power spectrum 
in a measurement of the redshift-space distortions, because their contamination could give 
rise to a systematic error when comparing with theoretical models. An assessment of the 
systematic error is also the purpose of the present paper. 

This paper is organized as follows. In section 2, we show the multipole power spectrum 
of the satellite LRGs and their contribution to the total LRG sample. In section 3, we 
introduce a simple halo model for a system consisting of central galaxies and satellite 
galaxies, then we show that the halo model with the LRG HOD explains the multipole 
spectrum. In section 4, we also discuss about the impact of satellite galaxies assuming 
a future survey targeting H-alpha emitters at high redshifts. Section 5 is devoted to 
summary and conclusions. Throughout the present paper, we adopt the Hubble constant 
Hq = 100/ikm/s/Mpc with h = 0.7 unless otherwise stated. 

2 Results of the LRG samples 

In this section, we demonstrate the contribution of the satellite LRGs to the multipole 
power spectrum. We here use the halo sample described in [11] using observed SDSS 
DR7-Full LRG sample in Northern sky (publicly available catalog prepared by [20]). The 
sample consists of 96762 LRGs with the magnitude —23.2 < Mg < —21.2 in the redshift 
range 0.16 < z < 0.47 (the mean redshift is 0.32) covering 1.44(Gpc/h)^ comoving volume. 
Halo is identified with the counts- in-cylinders techniques developed by [8]: two galaxies are 
considered neighbors when the transverse separation Ar_|_ < 0.8Mpc/h and the redshift 
difference Az/{1 + z) < 0.006 corresponding to the velocity difference 6vp = 1800km/s. 
The total number of halos is 92046. When the missing galaxies due to fiber collisions are 
taken into account, the total number of LRGs become 98991. If all of them are hosted 
by the same halos of the observed LRGs, the actual number of satellite LRGs becomes 
6945 (7%). Most of halos (95.5%) occupy single LRG (hereafter we call them "single LRG 
systems") and the rest of them contain multiple LRGs ("multiple LRG systems"). The 
multiplicity distribution of LRGs in halos is listed in Table 1 of [11]. For the multiple LRG 
systems, we choose the brightest LRG (BLRG) in each group as the central LRG and the 
rest of them are the non-brightest LRGs (NBLRGs), which we regard as satellite LRGs. 
Strictly speaking, BLRGs are not always central LRGs, and our satellite sample contains 
central LRGs to some extent. We have used different samples described in Table 1 to see 
the impact of satellite galaxies on the redshift-space power spectrum. 
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Name of LRG sample 


Total Number of LRGs 


Number of non-brightest LRGs 


All 


96762 


4716 


BLRG 


92046 





Single 


87889 





NBLRG 





4716 



Table 1: Properties of LRG samples: "All" include all of LRGs in the SDSS DR7 LRG 
Sample in Northern sky; "BLRG" includes the brightest LRG in each group and the 
fainter LRGs are excluded; "Single" includes LRGs in single systems only, and any LRGs 
in multiple LRG systems are excluded; "NBLRG" consists of LRGs in multiple LRG 
systems except for the brightest LRGs. 

Figure 1 shows the histogram of the number density of the galaxy samples as a function 
of the redshift z. 

2.1 Multipole power spectrum 

We adopt the method to measure the multipole power spectrum developed in [21]. For 
simplicity, we adopt the weight factor '0 = 1, and the parameter a = 0.1 for the random 
catalog (see [21] for details). The method doesn't take the window effect of the survey 
region into account, but it is demonstrated that the window effect in our method is negli- 
gible by comparing with the other method incorporating it explicitly [19]. We perform the 
multipole power spectrum analysis for each sample, whose results are shown in Figures 2 
and 3. 

Figure 2 compares the multipole power spectrum of the All LRG sample (black curve) 
and that of the NBLRGs (red diamond with large error bars). The squares with small 
error bars are the results in a previous work in [6], which are obtained from the LRG 
sample with 7150 square degrees sky coverage with the total number 100157 in the range 
of redshift 0.16 < z < 0.47. Thus the previous sample is almost same as the "All" LRG 
sample in the present paper. This figure shows that the amplitude of the correlation of 
the NBLRGs is quite large compared with the dominant component. 

Figure 3 shows the effect of the contamination of the NBLRG on the multipole power 
spectrum. The black curve is the results with the All LRG sample, the green curve is the 
one with the Single LRG sample, and the blue curve is the one with the BLRG sample. The 
difference between the green and blue is small, which means that the difference between 
the Single LRG sample and the BLRG sample is small. But the difference between the 
black curve and the blue curve is significant, which means that the contribution from the 
NBLRG sample is crucial though the fraction of the NBLRGs are small. This feature is 
significant for P2{k^ and P4,{k), especially. Thus, the contamination of the satellite galaxy 
is quite important in these multipole power spectra. 
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Figure 1: Histogram of the number density of the LRG samples. The black, blue, green and 
red curves show the number density of the All, BLRG, Single, and NBLRG, respectively. 

2.2 Impact on parameter estimation 

Here let us demonstrate the impact of the contamination from the satellite galaxies 
(NBLRGs) in an estimation of cosmological parameters. For simplicity, let us consider 
the simple model of the anisotropic power spectrum 

P(fc, /x) = (bik) + ffiYP^'-imkfia,], (3) 

where P^^{k) denotes a nonlinear matter power spectrum, Vlkfia^] is the damping factor 
due to the FoG effect and is the velocity dispersion parameter, for which we adopt 
the function V[k^av] = 1/(1 + Here we determined the bias b{k) so that 

the observational and the theoretical monopole spectra match. Then computed the chi- 
squared using the quadrupole spectrum by = ^i[P2^'^' (ki) — P2^''°'{ki)]'^/[AP2^''-{ki)]'^, 
where P2^^'{ki) and AP2^^'(ki), are the observed values and errors, and P^^^'^iki) is the 
corresponding theoretical value. See Ref. [5] for details. 

Figure 4 shows the 1 sigma and 2 sigma contours of Ax^ on the parameter plane 
and 7, where the growth factor and the growth rate are parameterized as 



Di{a) = a exp 



Jo « 



(4) 



dlogPija) 

fW = — Tj = ^m{ay, (5) 

a log a 

where ^^(a) is the matter density parameter at the scale factor a. Here we fixed the 
other parameters Ug = 0.97, flm = 0.28, fib = 0.046, as = 0.8 and assumed the cold dark 
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Figure 2: Multipole power spectrum Po{k), P2{k), P^^ik), and PQ{k) for the All LRG sample 
(black curve) and for the NBLRG (red diamond with large error bars). The squares with 
the small error bars show the results with the sample in a previous paper for comparison. 
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Figure 3: Same as Fig. 2 but comparing the results with the All LRG sample (black curve), 
the BLRG sample (blue curve), and the Single LRG sample (green curve). The red curve 
is obtained by summing each component of the BLRG sample, the NBLRG sample and 
the cross correlation component. The agreement between the red curve and the black 
curve shows a consistency of the computation. 
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Figure 4: Contour of Ax^ on a^, and 7 plane. The solid (dotted) curves are the 1 sigma and 
the 2 sigma contours with the power spectrum with the brightest (single) LRG sample, 
while the dashed curve is the same but with the All LRG sample. The left panel used 
the data in the range of wavenumbers 0.01/iMpc~^ < A; < 0.2/2,Mpc~^, but the right panel 
used the data in the range O.Ol/iMpc"^ < A; < 0.3/iMpc"\ 

matter model with a cosmological constant as the background universe model. In each 
panel, the dotted curve, solid curve, and the dashed curve are the Single, Brightest, and 
All LRG sample, respectively. The left (right) panel used the data with k < 0.2/iMpc~^ 
{k < 0.3/iMpc~^). The value 7 = 0.55 is the prediction of the model on the basis of the 
general relativity [22]. Though our theoretical model is very simple, the results clearly 
show that the contamination of the satellite galaxies (NBLRGs) significantly biases the 
parameter estimation. This figure also indicates that the results are influenced by including 
the brightest LRGs consisting of the multiple systems. 

3 Analytic model 

In this section, we consider the FoG effect of satellite galaxies based on the halo model 
picture [23, 24, 25]. In the halo model, the power spectrum of LRGs are decomposed into 
1-halo and 2-halo terms. Then we write the anisotropic power spectrum in the redshift- 
space consisting of the 1-halo and 2-halo terms, 

PLRG(fc,/i) = P'\k,fx)+P^\k,fx). (6) 

We here consider the sample which consists of the central galaxies and the satellite galaxies, 
and adopt the following expressions (7) and (13) for P^^{k,n) and P'^^{k,fi), respectively. 
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Figure 5: HOD for LRGs based on [8]. 



A brief summary of the derivation for a general case is described in the appendix (See also 
below for details). 

One-halo term is given by 



■n? 



dM 



dn 



dM L 



2(Ar,en)(iVsat)^cs(A;, ^l] M) + (iVsat(iVsat " l))Pss{K f^] M) ,(7) 



where we adopt the halo mass function dn/dM given by [26] and n is the mean number 
density of LRGs given by n = / dM{dn/dM)NnoD{M) and NnoD{M) is the halo occupa- 
tion distribution (i.e., the average number of galaxies inside the halo with mass M). We 
use the following form of the HOD of central LRGs and satellite LRGs [32] 



iVHOD(M) = (iV,en)(l + (iVsat)), 



(iVccn) 
(iVsat) 



1 + erf 



M -Mr 



log,o(M)-logio(M„ 



Ml 



(9) 

(10) 



where erf(x) is the error function. We adopt Mmin = 5.7 x IO^^Mq/H, a\ogM = 0.7, 
Mcnt = 3.5 X lO^^Mo/Zi, Ml = 3.5 x IO^^Mq/H, and a = 1 to match the HOD of SDSS DR7 
LRG catalog [8] as shown in Figure 5. Assuming the number of groups with A'sat satellites 
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is Poisson distributed [27], the averaged satellite-satellite pair number (A'satlA^sat — 1)) per 

halo goes to (iVcen)(iVsat)^- 

Central LRGs locate near the halo center and thus their velocity difference relative to 
the host halo should be small. Note that it is difficult to verify that each central LRG 
is located at the center of each halo in observational data. However, 20-40% of brightest 
LRGs are found to be off-centered (satellite) galaxies using lensing and cross-correlation 
analysis [11]. Therefore, large part of the NBLRGs are off-centered and their velocity 
should be the main source of the FoG effect. The functions pcs{k, /i; M) and pss{k, /i; M) are 
the Fourier transform of central-satellite and satellite-satellite distribution inside the halo 
with the mass of M, and the internal motion of satellite LRGs elongate the distributions 
in the line-of-sight direction. We assume that the internal velocity of the satellite LRGs 
has a Gaussian distribution determined by virial velocity as ai,^os{M) = (GM/2i?vir)^''^, 
in which the virial radius of the halo with mass of M is i?vir = (3M/47rpni(-z)Avir(-z))^^^ 
with Avir = 265 at z = 0.32. When the satellite motion is uncorrelated with each other, 
pcs{k, fi, M) and Pss{k, fi, M) are given by 



Pcs{k, /i, M) = unfw(^; M) exp 



:iii 



PUk,fi,M)=pi{k,fi,M). (12) 

We assume that the distribution of the satellite galaxies follows the NFW profile [28] and 
■Unfw(^) denotes the Fourier transform of truncated NFW profile [29]. 
The 2-halo term is given by 



d/n 

- / c^M— (6(M) + //i2)(Ar,,,)(l + (Ar3,,)p,3(A:,/i;M))M™i(A;;M) 
n / dM 



P^k), 
(13) 



where P^^{k) is the real-space non-linear matter power spectrum. Here we add the volume 
exclusion effect of halos u^oiik] M) in addition to the satellite distribution in order to 
include that two different halos cannot approach each other closer than a halo size. We use 
a Gaussian form Umi{k; M) = exp(— (aA;i?vir(^))^/2) and we choose the width parameter 
a = 2 to fit the observed power spectrum. 

Again we consider only the velocity distribution of satellite LRGs. We simply use 
the linear Kaiser formula [30] given by the term of {b{M) + //i^) with the growth rate 
/ = d\nD/dlna and the linear halo bias for the halo bias b{M) [31, 29, 25]. Without the 
FoG effect, that is ay^os = 0, we have P"^^ ~ (6efr + //i^)^-PNL and P^^ ~ iVih, where bes 
is the effective bias of LRGs given by b^s = f dM{dn/dM)b{M)NuoD{M)/n and A^ih is 
defined by Ni^, = J dM{dn/dM) [2{N^^^){N,J + (iVsat(iVsat - 1))] /n^- In this case, Nu, is 
a constant and we have Pj^ik) = for i > 2. 

Figure 6 compares the halo model predictions of multipole power spectra for All LRG 
sample with the observations. The model predictions well agree with the observations 
although we simply adopted the linear Kaiser redshift distortion and the linear halo bias. 
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Figure 6: Halo model prediction for the multipole power spectra ^0(^)5 P2{k), P^ik), and 
P(,{k) for the All LRG sample. In each panel, the dotted curve and the dashed curve are 
the 1-halo term and the 2-halo term, respectively, and the solid curve is their combination. 
The black circles are the observational data of the All LRG sample in Fig. 3. 
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Figure 7: Halo model prediction for the multipole power spectra Po{k), P2{k), P^^ik), and 
Pel^) for the NBLRG sample. In each panel, the dotted curve and the dashed curve are 
the 1-halo term and the 2-halo term, respectively, and the solid curve is their combination. 
Here the BLRG sample is assumed to be consisting of the central galaxies (35%) and the 
satellite galaxies (65%). The black circles are the observational data of the NBLRG sample 
in Fig. 2. 
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Figure 8: Halo model prediction for the multipole power spectra Po{k), P2{k), P^ik), 
and Pe{k) for the BLRG sample, which is written with just 2-halo term. Here the solid 
curve adopted the fraction of the central galaxies gcen'^*^'' = 1, while the dotted curve did 
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The satellite FoG effect in 1-halo term becomes significantly important at larger k and 
dominant in the multipole spectra for i > 4 even though the satellite fraction is only 
5%. The 1-halo term contribution causes a systematic bias in the measurement of the 
growth rate as shown in Figure 4 because the FoG effects from the 1 and 2-halo terms 
have different feature and the simple form of equation (3) cannot describe both of the FoG 
effects very well. 

The behavior of higher-order multipole spectrum is sensitive to the satellite FoG effect 
in one-halo term. In other words, the higher multipole spectra may be a good probe of the 
satellite fraction and the satellite velocity distribution. The one halo term contribution 
to the multipole power spectrum of (7) can be written as follows 

^^"'(^) = ^ / "^^dM [2(A^cen)(iVsat)g£(g) + (iVsat(iVsat " l))Qe{V2q)\ , (14) 



where we defined 



Qi{q) = ll'^d^ie^'^''^'C,{fi) 
and q = ay^os{M)k/ ^/2aH{z). Specifically, we have 



(15) 



Qo(g) = ^erf(g), (16) 



^{3- 2q^) 
4g2 ^ 8g3 

5(21 + 2g2) 2 30F(35 - 20g2 + 4g^^ 



Q2(.) = -7^e- +^^^^^^^-^erf(g), (17) 



QM = -^^^e- + ^''^ 

21(165 + 20g2 + 4g^) , 5^/^(693 - 378g^ + 84^^^ - Sg^) , 
Qd<l) = e + erf(g). (19) 

The error function has the asymptotic form erf(g) — )■ 1 for q ^ 1. Therefore, Qt{q) is in 
proportion to (— l)^/^g~^ in the limit g 3> 1, which explains the asymptotic behavior of the 
multipole power spectrum at the large wave numbers. The central-satellite contribution 
is dominant in the 1-halo term of LRG samples and thus kP}^ ~ '2kQ(,{q)fsa,t/n where /sat 
is the satellite fraction. When we use the values of cx^^ofF = 670km/s, aH{z) = 88/ikm/s at 
z = 0.32, /sat = 0.07 and n ^ 10~^(Mpc/h)~^ for the LRG sample, the large-scale limit of 
kP^^ goes to 230 for i = 0, 120 for £ = 2, 85 for £ = 4, and 74 for i = 6. Figure 6 shows 
that 1-halo term contribution approaches these values. 

Figure 7 shows the halo model predictions of the multipole power spectra for the 
"NBLRG" sample compared with the observed spectra. As 40% of BLRGs in multiple 
LRG systems are satellites (the number of multiple LRG systems AAmui is 4157), the same 
number of central galaxies are mixed in the NBLRG sample. Here we consider that 35% 
(=0.4jVinui/A/'sat) of the NBLRG sample are central galaxies and the rest of them are 
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satellites. Based on this assumption we write the one-halo and two-halo terms of the 
power spectrum of the NBLRG sample as 

jNBLRGn. ,,^ _ Dlh.NBLRGj^^ , D2h,NBLRG/ 

1 f ^ dn 

'sat 

'Hln^ (1 - Qt:^)Pcs{k, /i; M) + (1 - qt^^fPUk, /i; M)] , (21) 
dn 



P™LRG(j._^) _ pl'>™I-'»°(/fc,f,) + p2'>-™I'K°(/t,fl), (20) 



52h,NBLRG 



(/c,/i) 



rfM— (6(M) + //)(iVeen)(iVsat) 



2 



X {qt:^ + (1 - g^rn M))uUk; M)] (A^), (22) 



where gccn*^ is the fraction of central galaxies in the NBLRG sample. Here we set gcen*^ = 
Min(0.35, 1/ {Nsat)) so that qiln\Ns.at) does not exceed unity. The one-halo and two-halo 
terms for NBLRGs are plotted with the dotted curve and the dashed curve, respectively. 
The deviation at large scales [k < 0.2/i/Mpc) might come from the difficulty to make 
random samples because the number density of NBLRGs are very small. Nevertheless, 
the halo model qualitatively explains the observed multipole spectra at smaller scales. 
One-halo term (dotted curve) is a dominant contribution to the multipole spectra and 
reaches kP^^^^^{k) ~ 0(10^). 

The brightest LRG power spectrum does not have one-halo term because each halo 
contains one LRG at most. Then the BLRG power spectrum is written only with the 
two-halo term. This can be clearly seen in that the BLRG multipole power spectra with 
i = 4 and 6 are consistent with zero, which indicates that the 1-halo term contribution 
becomes negligibly small. However, several observations indicate that some fraction of 
BLRGs are satellite or off-centered galaxies (e.g., [35]), which causes the FoG effect [10]. 
When the fraction that BLRGs locate on the mass center of their host halos is gcrai"^*^'', 
the BLRG power spectrum is given 



1 /" lin 

- dM —{b{M) + ff,'){N^^^) 



(e^^^^ + (1 - e'^^^);^cs(fc, /i; M))S.o,] ' P^\k). (23) 

As shown in [11], the lensing and cross-correlation measurements indicate 20% fraction 
of BLRGs are off-centered. Figure 8 shows the model predictions of the BLRG power 
spectra with qcen^'^^ = 1 (all of BLRGs are centrals) and gcen"^*^^ = 0.8 (20% of BLRGs 
are satellite). The figure shows that their agreement is better when the satellite FoG 
effect is included, which indicates that even BLRG sample has significant FoG effect on 
the multipole power spectra. Note that the result may change if we take into account the 
nonlinearity in the galaxy biasing. The linear Kaiser formula is the simplest model, then 
more careful analysis will be necessary using the sophisticated perturbation theories as 
well as numerical simulations ([36, 37, 38, 39]), though such analysis is beyond the scope 
of the present paper. 
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4 Forecast on multipole power spectra for Ha emit- 
ters 



Main targets of high-redshift (z = 1 ~ 2) galaxy surveys, planed in such as Subaru/PFS 
[40] and Euclid [41], are Ha emitters (HAE). The relation of HAEs to halos are less known 
observationally, and will be more complicated than that of LRGs. We use the following 
form of HOD based on the sample of 370 HAEs at z=2.23 detected in Hi-Z Emission Line 
Surveys (HiZELs) [42] 



(^cfn ) = m - exp 



1 + erf 



(logio(M) - logio(M^i,,,))- 
logio(M) -logio(M^in,c 



(iVsat) = 



1 + erf 



C" log M.c 
logio(M) - log^o(Mmm,s 
"^log M,s 



M - M„ 



(24) 
(25) 



Here the central HAE distribution is described with Gaussian and smoothed step-like 
components with their amplitudes determined by the normalization factors Fa and Fi,. 
The typical mass and the dispersion are parametrized with Mmin^c and o"iogA//x, respectively. 
Satellite HOD is described with a smoothed step-like component multiplied by power-law 
with scaling of a, the typical satellite mass Mjnin,s and the amplitude Fg. The values of 
HOD parameters for different luminosity samples are listed in Table 2. The plots of HODs 
are shown in Figure 9. 



Luminosity 


L > lO^^erg/s 


L > lO^^erg/s 


Fa 


0.4 


0.0035 


Fb 


0.33 


0.06 


Fs 


0.1 


0.001 


M^i„,,(10i2MQ//i) 


0.12 


1 


log M.c 


0.14 


0.22 


M^i„,3(10i2MQ//i) 


0.5 


1 


Clog M,s 


0.24 


0.24 


a 


1 


0.8 



Table 2: HOD parameters for HAEs 

Figure 10 shows the comparison of multipole power spectra for central HAEs with those 
for all HAEs including satellites at z = 2. The FoG effect on the HAE power spectrum from 
satellite is much smaller than that on LRGs power spectrum because the typical halo mass 
of HAEs is much smaller than that of LRGs. In our halo model, averaged virial velocity 
of halos hosting HAEs is 170km/s, while those hosting LRGs are 670km/s. For the faint 
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Figure 9: HOD for Hq, emitters based on [42]. 



HAE sample, the contamination of the sateUite changes the higher multipole spectrum at 
a few percent or 10 percent level depending on the wave number, while the effect becomes 
smaller for the luminous HAE sample because the satellite fraction decreases. FoG effect 
for HAEs are expected to be much smaller than LRGs, however, upcoming galaxy surveys 
are expected to measure the growth rate measurement at the percent-level accuracy and 
thus it is still important to estimate the systematic errors of the FoG effect. 



5 Conclusion 

In the present paper, we have investigated the influence of the satellite galaxies on a 
measurement of the redshift-space distortion. We have found the following points for 
the first time. First, the satellite galaxies significantly contribute to the higher multipole 
power spectrum though the fraction is small. Second, the contribution of the satellite 
galaxies to the higher multipole power spectrum is explained by a simple halo model, and 
the one halo term in the halo model is a dominant contribution. These findings are based 
on the SDSS LRG sample, and the contribution from satellite galaxies may depend on 
the HOD of a galaxy sample. This means that an uncertainty of the HOD might give 
rise to a systematic error in measuring redshift-space distortion when satellite galaxies 
are contaminated. For the Ha emitters, which are the target galaxies of the PES redshift 
survey and the Euclid redshift survey, we have shown that the satellite's contribution to 
the redshift-space distortion is much smaller than the case of LRGs, because the host 
halo mass is small. The results are based on the HOD of the Ha emitters at the redshift 
z > 2, it would be interesting to investigate how the results change depending on the 
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Figure 10: Multipole power spectra for Hq emitters with L > lO^^erg/s (left panels) and 
with L > lO'^'^erg/s (right panels). 
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redshift especially in the lower redshift regions. Combination with weak lensing survey 
might help to resolve the uncertainty in HOD [44, 45, 46]. However, we should investigate 
the systematic error from the HOD uncertainty in the cosmological parameter estimation 
at a quantitative level. These problems are left as future work. 
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A Derivation of power spectrum 

In this appendix, we derive a general expression of the multipole power spectrum in the 
halo model, which gives the grounds to adopt the expressions in section 3. Following the 
halo model approach, the correlation function is written as the sum of the 1-halo term and 
the 2-halo term. The power spectrum is the Fourier transform of the correlation function, 
then the power spectrum is also written as the combination of the 1-halo term and the 
2-halo term. We start with the real-space power spectrum in a halo model presented in 
Ref. [33], 



where we defined 



dM^^{N,en) 



P^'\k) = - 



dM 



dM 
dn{M) 



2{Nsat)u{k\M) + {NUNsat - l))u{k\Mf 



dM 



{N,,n){l + {Nsat)u{k\M))h{M) 



Pm{k), 



(26) 



(27) 



(28) 



and u{k\M) is the Fourier transform of the density profile of galaxy distribution. We 
assume that the galaxy density profile is the same as the dark matter density profile. For 
the NFW density profile, we have [25] 



u{k\M) 



ir<7\ 



3 
s 



M 



d^xp{x\M) 
sm{krs) [Si{[l 



(29) 



c]krs) - Si{krs)] 



sin ckr^ 



1 + c)kr 

+ cos(fcrs) [Ci{[l + c]krs) - Cz(A;r,)] 



(30) 
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where 



CAx) = - 



cos t 



S^{X) 



sint 



dt. 



(31) 



The redshift-space power spectrum of the halo model may be evaluated as follows. 
Tinker investigated the formulation for the redshift-space correlation function in a halo 
model [34], in which the redshift-space correlation function is obtained by 



e{r)V{v,)dv, 



(32) 



where ^^{r) is the real-space correlation function, r± is the projected separation, ry is 
the line of sight separation, = + z'^ and = H{r\\ — z), V{vz) is the probability 
distribution function of the galaxy pairwise velocity, and H is the Hubble parameter. 
Assuming the Gaussian distribution function, V{vz) = (7r^/^2(T,u)~^e~'^^/^°"", we may write 
the redshift-space power spectrum in the form 



where 



V? dM ^ ' 



2{N,at)u{k\M)e-''"^"^'"'^''"''" 



and 



n dM ^ ' 



X 1 + 



(iV,,i)M(A;|M)e^'^''^'^'/2'^'^') (6(M) + ///^ 
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